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Abstract 

We present here a compactly formulated application of the previously posted general formalism 
of the reflection of Gaussian beams at a dielectric interface (arXiv:0710.1643v2 [physics.optics]). 
Specifically, we calculate the Goos-Hanchen shift near Brewster incidence, for an air-glass plane 
interface. 
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I. INTRODUCTION 



When a beam of light impinges upon a plane interface separating two transparent media, 
it produces reflected and transmitted beams. In 1815 the Scottish physicist David Brewster 
discovered the total polarization of the reflected beam at the angle 9b since named after 
him jl|. From his observations he was also able to empirically determine the celebrated 
equation, known as Brewster's law, tan^^ = nxjn2^ where and are the respective 
refractive indices of the two media. 

In this work we calculate the Goos-Hanchen shift occurring near Brewster incidence at 
an air-glass plane interface, for an incident Gaussian beam. 

II. THEORY 

Consider a monochromatic beam of light incident upon a plane interface that separates 
two homogeneous and isotropic media. The first medium, say air, has refractive index n< 
and the second medium, say glass, has refractive index n>. With n = ^ we denote the 
ratio between the two refractive indices. Here n can be either a real or a complex number, in 
the latter case at least one of the two media exhibits absorption. Without lack of generality, 
we assume that the beam meets the interface coming from the air side. Thus, it will be 
convenient to take the axis z of the laboratory Cartesian frame K = {O, x, y, z) normal to 
the interface and directed from the air to the glass. Moreover, we choose the origin O in 
a manner that the plane interface has equation z = 0. The air-glass interface, the incident 
and the reflected beams are pictorially illustrated in Fig. 1. 

In addition to the laboratory frame, we use a Cartesian frame Ki = {0,Xi,yi, Zi) attached 
to the incident beam and another one Kr = {0,Xr,yr, Zr) attached to the reflected beam. 
Let ko = koZi and k denote the central and noncentral wave vectors of the incident beam, 
respectively, with |k| = |ko| = kg. We choose the laboratory frame K in such a way that 
Zj = xsin^^ + zcos^^. In this manuscript with either u or we denote a real unit vector 
directed along the Cartesian frame axis u, where u G {x, y, z}, and a G {i, r}. 
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FIG. 1: (Color online) Geometry of beam reflection at the air-medium interface. 9b is the Brewster 
angle. 

The electric field of the incident beam at the air side of the interface {z < 0), can be 
written in the angular spectrum representation ^ as 

E\r) = J2 [[ ^' ^)^A(f/, \/)e^(^^'+^^'+^^-) dUdV, (1) 

A=l 



-oo 



where = k • Xj/fco, V = \<. - ji/ko, and W = {1 — U'^ — V'^Y^'^. Moreover, we have defined 
Xa = koXa, Ya = koZa, Za = k^Za, with a G {«,'r}, and r so . The polarization unit basis 
vectors {eA}Ae{i,2} have been chosen as 

62 X k „ z X k 

|e2 X k| |z X k| 

where the symbol "x" denotes the standard vector product in . Here z is a real unit 
vector directed along the laboratory axis z and, by definition, 

k = fco {U±i + Vy, + W7.i) 



yikx + yky + 7.J kl - kl - kl. (3) 
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By substituting Eq. ([3]) into Eq. ([2]) we obtain 

. , . „, [UWc + {W^ + V^) s]±, + V {Wc - Us) y, - [{U^ + V^)c + UWs] z, 
ei(ty, V, 6) = — , {A) 

[V^ + {Uc + Wsf] ' 

\V^ + {Uc + Wsf] 

where we used the shorthand c = cos^ and s = sin 6', and 6 is the central angle of incidence 
defined as: 6 = arccos(zj ■ z). For well collimated paraxial beams {U"^ + V^^ <^ 1) the 
expressions above reduce to 

ei ~ Xi + 1/ cot e ~ZiU, (6) 
ea ~ -Xi V cot ^ + Yi - Zi V, (7) 

63 ~ Xi [/ + + Zi, (8) 

where 63 = k/fco = k. From Eq. it follows that ei lies in the plane containing both 
the wave vector k and z, usually denoted as the plane of incidence with respect to k, while 
62 is orthogonal to such a plane. Both ei and 62 are orthogonal to k by definition, and 
{ei,e2,k} form a complete, orthogonal basis in M.^. Conventionally, a plane wave whose 
electric field is parallel to either ei or 62, is referred to as either a TM or a TE wave, 
respectively. The symbols P for TM and S for TE, are also widely used. In Eq. ([1]) the 
functions Ex{U, V) determine the shape and the polarization of the beam. These amplitudes 
are always expressible as E\{U, V) = A{U, V) a\{U, V), where A{U, V) and ax{U, V) are the 
scalar and the vector spectral amplitudes of the field, respectively. The first determines the 
spatial characteristics of the beam, while the second sets the polarization of the beam j^. 
Here we consider a collimated, monochromatic beam, with a Gaussian spectral amplitude 

A{U, V) = exp ^ j exp (zWD) , (9) 

where ^0 = 2/(A;oWo) is the angular spread of the incident beam with a minimum spot size 

(waist) Wq located at = —D In order to determine the vector spectral amplitudes 

ax{U,V) of the incident beam, we assume that the beam has passed across a polarizer 

plate perpendicular to the central wave vector koZi of the beam. Let f = /pXi + fsYi 

denotes a complex-valued unit vector that represents the orientation of the polarizer, with 

l/pp + l/pp = 1. Then, we determine the amplitudes ax{U, V) by imposing 

2 

J] e,(f/, V, e)ax{U, = f - k(k ■ f) , (10) 

A=l 
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where we used the polarizer representation given in Ref. 4]- Since the completeness of the 
basis {ei, 62, k} implies, for any vector v, the validity of the following relation 

ei(ei ■ v) + 62(62 ■ v) = V - k(k ■ v), (11) 

then from Eq. ffTOl) it immediately follows that 

ax{U,V)=e^iU,V,e)-i. (12) 

Thus, from Eqs. (jUlS]) and f|T2|) we obtain 

W {Ufp + Vfs) cose -[UVfs-{W^ + V^) fp] sin 6 

[V^ + iUcose + WsinOf] 
{Ufs-Vfp)cose + Wfssm e 

[V^ + {Ucose + Wsmef] 
that reduce, for paraxial beams, to 



aiiU,V)c^fp + fsVcot9, (15) 
a2{U,V)c:^fs-fpV cote. (16) 

When the beam is reflected at the interface, each plane wave mode function 

exiU, V, ^)e*(t'^>+v'^»+w/z.) ^ eA(k)e^''■^ (17) 

changes according to 

where ri(k) and r2(k) are the Fresnel reflection amplitudes for TM and TE waves, respec- 
tively 5I, 

n^kz + kl kz + kl 



and k = k — 2 z (z ■ k) is set by the law of specular reflection [g], while the unit vectors eA(k) 
are defined as in Eq. ([2]) with k i-^ k. In Eqs. ( 1191) k\ is the z-component of the wave vector 
inside the glass, namely 

k\={n^kl-kl-k^"\ (20) 
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It is worth noting that here ka:,ky,kz are the Cartesian components of the wave vector k 
with respect to the laboratory frame K, while in Eq. ([1]) the integrations are performed 
with respect to the variables U and V which are the transverse Cartesian components of the 
wave vector k with respect to the incident-beam frame Ki. Therefore, it will be useful to 
express ri(k) and r2(k) in terms of U and V. From Eq. ([2]) it straightforwardly follows that 

rr I vv cos o — u sni at— irr — v ~ — (vv sni a -\- u cos (7)^] ^ 

^^(u, V) = ^ y^, (21) 



r2iU,V) 



{W cosO-U sin 6) - [n^ -V^ -{WsmO + U cos Of] 
n2 {W cosO-U sin 9) + [n^ -V^ - {Wsine + U cos Of] ^''^ 
cos e - f/ sin e - - _ gin ^ + f/ cos Qf\ ^'"^ 



(22) 



Wcoi^e-Uis,ine+[n^ -V^ -{Wsixie + U cos dy] 
It is easy to check that ri(f/, V) and r2(f/, V) reduce to the ordinary Fresnel coefficients for 
[/ = and 1/ = 0: 

,2 Ml/2 



ri'(^) = ri(0,0) 



71^ COS 



71^ COS 



e - -sin^^) 



sm' 



1/2 ' 



cos — (ra^ — sin^ 6') ' 



cos 9 + (n^ 



(23) 



(24) 



sin^^)''' 

In the remaining of this manuscript, we shall often benefit from the following relations 
satisfied by the refiection coefficients defined above: 



drx 



dU 



and 



u=o,v=o 



d'^rx 



dri 



drx 



dV 



(7=0, y=o 



u=o, v=o 



0, 



d'^rx 



dUdV 



0, 



(25) 



(26) 



(7=0, y=o 



[7=0, y=o 
where A G {1,2}. 

From Eq. (fT8|) it follows that, after refiection, the electric field of the beam can be written 



as 



A=l 



E^(r) ^ E^(r) = J2 ^^^i^^^i^, V)Ex{U, yy'^-v^r^-VYr^-wz.) ^27) 



where eA(k) = Qx[—U ,V — &) , namely: 

^WC + (iy2 + -1/2) s\%-V {y^C - Us) Yr + [(f/2 + V^)C+ UW s] Zr 



ei(k) 
62 (k) 



[V^ + iUc + Wsf]'^'' 

Vc^ir + {UC + Ws) yr - VS Zr 

[V^ + {Uc + WsfY^^ 



(28) 
(29) 
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where we used again the shorthand c = cos^ and s = sin 6'. In Eq. (127|) we have exploited 
the fact that by definition 

k-r = Kx + kyy~ k,z = -UXr + VYr + WZ^, (30) 

where the latter equality is written in terms of the Cartesian coordinates of the position 
vector r with respect to the refiected beam reference frame K^. 

If the air-glass interface would behave as and ideal refiecting surface characterized by 
wave vector- independent refiection amplitudes ri(k) = 1 and r2(k) = — 1, then the refiected 
beam were just the mirror-image of the incident one [7j]. However, in the real world, as 
a result of the polarization and wave vector dependence of the Fresnel amplitudes rA(k), 
non-specular refiection phenomena occur, the most prominent of which are the so-called 
Goos-Hanchen (GH) jsl and Imbert-Fedorov (IF) shifts that amount, respectively, to a 
longitudinal and a transverse displacement of the refiected beam with respect to the mirror- 
image of the incident one. Such displacements can be assessed by measuring the position of 
the center of the refiected beam with a quadrant detector centered at a;,. = 0, = along 
the reference axis attached to the refiected central wave vector ko = /co^r- A quadrant 
detector has four sensitive areas each delivering a photocurrent when illuminated. The 
difference between these photocurrents is proportional to the displacement of the barycenter 
of the beam intensity liX^-^Y^^Z^^ relative to center of the detector. In other words, this 
displacement is proportional to the first order moments (X) = (Xr)xr + (^r)yr of the 



intensity distribution function of the beam jlO| : 



(X) = . (31) 

y ^ 1[X,., Yr, Zr)AX,AYr 
— oo 

In order to evaluate (X) we need to know the intensity /(X^, yj,, Z^) that, apart from an 
irrelevant proportionality factor, can be defined as 

/(x^, y^, ^j.) = |E^(Xr, y^, ^^) 1^. (32) 

Thus, we must calculate the double integral in Eq. (l27j) . To this end, we exploit the fact 
that for a well coUimated beam 6q <^ 1, and that Eq. ([9]) implies that A{U, l^) ~ outside 



the paraxial domain V = {U,V : U'^ + V"^ 1}. In this domain 

W={l-U'- V'Y^' ^ 1 - (33) 
and we can rewrite Eq. fl27|) as 

oo 

E^(r) = e^(^^+^) [feiU, v^)e-^[A+^(^.+^)le^(-c/^^+v^>'0 d[/ dV, (34) 



-oo 



where we have defined A = 2/6q = koL, with L equal to the Raleigh range of the beam j^. 
Equation is still exact, and it defines £{U,V) as 

eiU^V) = X:e.(k)r,(f/, V)a.(f/, V)e^("^^"^K^-"^)], (35) 

A=l 

which can be evaluated within the paraxial domain V via a Taylor expansion of the form 
£{U, V) ~ £(0, 0) + {U£u + V£v) + ^ {U^£uu + 2UV£uv + V^£vv) +■■■, (36) 

where we used the obvious notation £u = d£{U,V)/dU\jj^Qy^Q, and so on. Usually, to 
calculate both GH and IF shifts, first order Taylor expansions is enough. However, as we 
shall see soon, at Brewster incidence it becomes necessary to keep second order terms to avoid 
divergences in the expressions of the shifts. Substitution of Eq. (!36|) in Eq. permits 
the analytical evaluation of the Gaussian integrals; this leads to the following expression for 
the electric field of the refiected beam: 

E^(r) ^ ^(r) {±,.E^^ + y,.El + z.EfJ , (37) 

where 



1 fiX'^ + Y'^ 
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is the scalar amplitude of a fundamental Gaussian beam, and 



X2(4-2rp) Y^[r'p-2{rp + rs) cot e]coie 
+ /p r./^ — nrT^5 — ^ 



2{Z-iKf 2(Z-zA)2 
+ 2r5 - 2 (rp + r^) csc^ 6* + cot 9 
2{Z -iA) 



p- 



„ XY [2rp + rs + rs cos(20) - {r'p + r^) sin(2^)] csc^ 9 
+ 2(Z-^A)2 ' 



Z - iA Z -ik 

„ XV^ „ [2rp - 2 (rp + rs) csc^ 9 + r'g cot 0] 



2(Z-zA)2 2(Z-iA)2 
+ 2rp — 2 (rp + rs) csc^ ^ + r^ cot 9 



s- 



2{Z -iA) 



XY [2rp + rs + rs cos(2^) - {r'p + r'g) sm{29)] csc^ 9 
"^^ 2(Z-zA)2 ' 

_ . f Yrs 

— ~ Jp-^ — ~ js- 



Z-iA '^Z-iA 
Y^rp + rs) cot 9 X\'p XY [r's - {rp + rs) cot 9] 

{Z~iK)^ ^^{Z~iA)^^^ {Z-iAf 
_^r/p + irp + rs)cot9 

Z-iA ^ ^ 

For sake of clarity, in the formulas above we have omitted the subscript "r" from the coor- 
dinates X, Y, Z, and we have used the shorthand 

r^:=r°(^) r'^--=^{(^) r'i:=^{9), {Ae{P,S}). (42) 

Moreover, as the variable D appears always in the form Zr + D, in the equations above with 
Z we denoted Zr + D, which amounts to a trivial re-definition of the origin of Kr. 

It is easy to see that the expressions for the electric field obtained above take explicitly 
the form of a power series expansion in the parameter 9q if we redefine the coordinates as 

2 2 2 

X = koXr = —C, Y = koVr = -^V, Z = ko{zr + D) = -^C, (43) 

where ^ = Xr/wo, rj = Vr/wo, and ( = {zr + D/ko)/L. After this rescaling, Eq. fl37j) takes 
the form of a power series: 

E^(r) ^ Eo^(r) + ^oEf (r) + 9lE^{r), (44) 
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where we have omitted an irrelevant overall multiplicative factor 9q. Finally, from Eq. 
the field intensity may be straightforwardly calculated as 

/(r) ^ |Eo^|2 + 9, (E^ ■ Ef * + c.c.) + Ul (|Ef 1^ + 2^^ ■ Ef + c.c.) + O{eo)\ (45) 

where "c.c." stands for complex conjugate. The explicit expression for /(r) is quite cumber- 
some and it will not be reported here. 

At this point, we have all the ingredients to calculate Eq. ( 13T|) that gives 



z 


l/p 


Vpr^ + 1/5 


1 1 
rs rg 


A|/p 




^ + ep) + \fs 





(^r) = - T^^^^^^T^^H^^k^T^^ (46) 



where ep represents the contribution of second order terms in the Taylor expansion and it 
is defined by 

ep = ^ [r'p + rprp-rl + rp r'p cot 6 + [r], - r|) csc^ 6] . (47) 

Here es is obtained from ep by interchanging the indices P and S. Note that for a TM- 
polarized beam at Brewster incidence rp = and fs = 0, and the denominator of Eq. (H7|) 
remains non zero only thanks to ep. 

Equation (H6ll shows that the distance from the beam center to the reference axis Zr grows 
linearly with Z as (Xr) = ZQ, thus defining unambiguously the angular shift of the beam 
equal to Q = d {X^) jdZ This definition is purely analytical and therefore, contrarily 
to the geometric one adopted by several authors 12|, llSj, it is always valid, even in the case 
of strong deformation or splitting of the reflected beam. 

In our experimental setup, beam reflection occurs at the front surface of a BK7 prism 
with refractive index n = 1.51031 at 826 nm, which corresponds to a Brewster angle Ob = 
arctann = 56.491°. For a TM-polarized incident beam (/p = 1 and fs = 0), Eq. (H^ 
becomes 



e = (48) 



which shows that Q = at Op where rp = and ep 7^ 0. However, since ep oc 1/A, and 
rp (X —{9 — Op) nearby 9p, then in this region Eq. (HHil reduces approximately to 

Aie-9py + a' ^ ^ 
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where 



1 



a 



+ 



2 ■ (l + n2)^' 

Since a does not depend on A, it is easy to see from Eq. 
X = 9 — 9bi then the following scaling property holds: 

v/A/(x/v^) = 



(50) 



£]) that if we put O = /(x), with 



(51) 



Thus, there exists an angle Om = Ob + a/o/A, {6m = 6b — A/a/A) close to where 9 
reaches a maximum (a minimum) approximately equal to l/(2-\/aA) [— l/(2\/aA)]. Since 
A = {koWo)'^/2 = 2/6*0, where 6*0 is the angular spread of the incident beam [2], then the 



maximum angular displacement occurring at 6m = 63 + 6Q\foij2 ~ 6*^ + 0.54 6*0 will amount 
to l/(2VaA) = ^o/VSa ~ 0.46 ^0 < ^o- Moreover, for \6 - 6b\ > ^0, Eq. (|49]) furnishes 

^0 



e 



6r<<'^ 



(52) 



2(t7 - OB) 

which is a signature of the sub-diffractive nature of the phenomenon. In Fig. 2 approximate 
expression (119!) is compared with the exact result (148!) for a beam waist Wq = 30 /im. The 
agreement between the two curves is very good and we have verified that it increases for 
increasing wq. 



wq = 30 fjm 



0.004F 




-0.004 



1 Rp dRp 
A Rp- + ep 



6 - e„ 



A (0 - BbT + a 



-1 1 

e - Ob (degrees) 



FIG. 2: Geometry of beam reflection at the air-medium interface. 
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III. ALTERNATIVE ROUTE 



Both Goos-Hanchen and Imbert-Fedorov shifts can be calculated in an alternative manner 
that displays the both spatial and angular characters of the shifts. The starting point is Eq. 
we rewrite here as 



E^fr) 



oo 

J I E{U, ^)e^^^e'(-t^^'-+^^'-) dUdV, (53) 



where Z = Zr + D and 

2 

E{U,V) = J2^x{k)r^{U,V)Ex{U,V). (54) 

A=l 

After a straightforward calculation, it is not difficult to prove the validity of the following 
formulas: 

oo oo 

jj\E^{X,Y,Z)\^ dXdY = Jj\E{U,V)\^ dUdV, (55) 



oo oo 
^R/ V V ri'W'i AV — „• / / TP* ATT A-i/ 'V II lIPlS 



X|E"(X,y,Z)|' dXdF = -i // —■E*dUdV-Z // —\E\'dUdV, (56) 



Y\E^{X,Y,Z)\^ dXdY = 1 II ^■E*dUdV + Z jj ^\E\^ dU dV. (57) 



Thus, we easily obtain 



[f^.E* dUdV ! ! ^\E? dUdV 



JJ dU JJ W' 
{X) = Z^ , (58) 

j j\E{U,V)\^ dU dV j j\E{U,V)\^ dU dV 



oo oo 

ill^-E*dUdV !!—\E\'^dUdV 



dV JJ W 
(Y) = ^ + . (59) 

jj\E{U,V)\^dUdV jJ\E{U,V)\^ dUdV 
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The equations above show clearly the spatial and the angular contributions to the shifts. The 
angular part is the part proportional to Z . It is interesting to note that the Z-dependence 
is strictly linear, as these equations are exact. Moreover, if we remember that U = k^^ and 
W = kz^, then it is obvious that 

^ = tan 6*2;, and ^ = tanOy. (60) 



From Eq. (MJh it immediately follows that 

d{X) 
dZ 

d{Y) 
dZ 



(tane,)~-(0,), (61) 



(ts^ney) (6.^) . (62) 
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